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The nucleation process of quark matter in both cold and hot dense hadronic matter is investigated 
using a chiral approach to describe the quark phase. We use the Nambu-Jona-Lasinio and the 
Chromo Dielectric models to describe the deconfined phase and the non-linear Walecka model for 
the hadronic one. The effect of hyperons on the transition phase between hadronic and quark matter 
is studied. The consequences of the nucleation process for neutron star physics are outlined. 
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I. INTRODUCTION 

In the last few years the physics of dense matter is having a growing interest due to its apphcations in both terrestrial 
(heavy ion collisions) and astrophysical (supernovae explosions and neutron stars) systems. In particular, the question 
of the existence of quark matter in the universe is still an open problem. Neutron stars are among the densest objects 
in the universe, and, therefore, they arc probably one of the best candidates to host a quark dcconfincd phase in their 
cores. 

Quark matter nucleation in neutron star has already been discussed in several previous works. Zero temperature 
nucleation was studied in Refs. [iHali while in Refs. J10l - [l5| thermal nucleation was considered. In all these papers 
the MIT bag model [16[ has been used for the description of quark phase. Recently, the nucleation rate and the 
nucleation time have been calculated at zero and finite temperature within the same quark model [iJ: UM- This 
calculation has allowed to follow the thermal evolution of neutron stars from the young and warm proto-neutron stars 
to the cold (T = 0) and dclcptonizcd neutron stars. The crossover temperature above which the thermal nucleation 
dominates quantum nucleation was calculated and the consequences of the possible quark nucleation for the physics 
and evolution of proto-neutron stars were discussed. It was shown that proto-hadronic stars with a gravitational mass 
below a critical value will survive the early stages of their evolution without decaying into a quark star. The effect of 
neutrino trapping on quark nucleation has been considered explicitly in Refs. |17| and |15| . 

In the present paper we will investigate how the nucleation process depends on the quark model chosen to describe 
the dense quark phase. There are several chiral quark models in the literature that successfully describe the low 
energy properties of mesons and hadrons. In particular, we will consider two quark models which contain explicitly 
the chiral symmetry: the Nambu-Jona-Lasinio (NJL) model [I^ and the Chromo Dielectric Model (CDM) [13. The 
NJL model (see [20| for a recent review) has been extensively used to study quark matter in /3-equilibrium and quark 
stars |20l - [25| . The parameters of this model are fixed by low energy scattering data of mesons. Although, the NJL 
model contains dynamical chiral symmetry breaking explicitly, an important symmetry of QCD, it is unable to explain 
confinement. This feature is maybe the main drawback of the model. 

The QCD phase diagram determined with the NJL model, shows that, at very high densities and low temperatures, 
the most favourable deconfined phase could be in a color superconducting state in which two flavors or all three quark 
flavors are paired |26l . |27[ . High densities and low temperatures are the conditions which are expected to occur at the 
core of a neutron star some minutes after its birth. In the present paper we do not take into account the possible 
formation of a color-superconducting state of dense quark matter and we use the standard version of the NJL model 
without superconductivity |28| . Recently the deconfinement phase transition in proto-neutron stars within the SU(3) 
NJL model including 2SC colour superconductivity has been discussed by the authors of Refs. [231 and |17| . 

The CDM has been used to study the static and dynamical properties of the nucleon, and describes the confinement 
of quarks through their interaction with a scalar field x which represents a multigluon state and produces a density 
dependent constituent mass. Quark matter has been analysed within this model in several papers i36. 43;. |46| . 

The present work is organized as follows: first we give a brief description of the approach used to model the hadronic 
phase; next we review the main features of the NJL and Chromo Dielectric models and present the results for quark 
matter nucleation in both cold and hot dense hadronic matter. Finally, the consequences of the deconfinement phase 
transition for neutron star physics are discussed. 

II. EQUATION OF STATE 

In this section we present a brief review of the models used to build the equation of state (EOS) of stellar matter. 
The hadronic matter is described within the non-linear Walecka model (NLWM) [43] and quark matter within the 
NJL and the CDM models. We take ?i = c = 1 in all the expressions. 

A. The non-linear Walecka model 

The Lagrangian density, including the baryonic octet, in terms of the scalar a, the vector-isoscalar w^ and the 
vector-isovector p^ meson fields reads (see e.g. [^2, 31, 33]) 
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and Mg = Mb — g^BO' is the baryon effective mass. The quantity tB designates the isospin of baryon B. The niesonie 
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where the sum is over electrons, muons and neutrinos for matter with trapped neutrinos. In uniform matter, we get 
for the baryon Fermi energy epB = guiBi^o + gpBhBPoz + ykp^ + M^"^. 

We will use the GMl parametrization of NLWM due to Glendenning and Moszkowski [30,[33]- The nucleon coupling 
constants are fitted to the bulk properties of nuclear matter. The inclusion of hyperons involves new couplings, which 
can be written in terms of the nucleonic ones: go-y = Xa- gcr, gu>Y = x^ g^j, gpv = Xp gp. In this model it is assumed 
that all the hyperons in the octet have the same coupling. Measured neutron star masses can be used to restrict the 
possible ranges of variability of the hyperon couplings J30ll33l|. In this work we will consider Xa = 0.6 and x^ = 0.7. In 
addition, following Ref. |30[ we will take Xp = Xa, whereas the binding energy of the A in symmetric nuclear matter, 
-Ba, 
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-28 MeV = x^ g^ wq - x^r gacr 



(10) 



is used to determine Xuj in terms of Xa- Notice that the case with Xa ~ 0.6 produces stars with a larger hyperon 
population (for a given stellar gravitational mass) with respect to the case x^ = 0.7 (7|, |33|. In addition to these two 
parametrizations for hyperonic matter, we will consider the case of pure nucleonic matter (hereafter called np matter). 
At zero temperature the equations of motion for the various fields can be derived using the Lagrangian of Eq. ([T]) 
and the well known Euler-Lagrange equations. At finite temperature we use the following grand canonical potential 
per unit volume: 
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where /? = 1/T, E* = (kf + Af*^)^/^ and the effective chemical potential of baryon i is given by 

z/j = ^ii - g^ujo - T3i gp po3 ■ (12) 

The sum is extended to all baryons of the baryonic octect. The equations of motion are then obtained minimizing fl 
with respect to the fields a, cjq and po3- The baryonic pressure and energy density are given by: 

Ps = -O , £ = -Pb + Y1 ^^"^ + ^bT (13) 

where rii is the particle number density of the baryon i and sb is the total baryonic entropy density: 

The total energy density and pressure are finally obtained adding the lepton contribution to Pg, eb and sb- 

B. The Nambu— Jona— Lasinio model 

The three flavor NJL model [l^l has been widely used by many authors to describe quark matter [20l - |25j . In this 
paper we adopt the Lagrangian density of Ref. [23 : 
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C = q{t-ff,d^ - mo)q + G^[(gAfeg)2 + [qij^Xkqf] ~ K[detf{q{l + 75)9) + detf{q{l - 75)9)] . (15) 

fc=0 

In the above expression, q denotes the three flavor quark field: q = {u,d,s) and toq = diag{mQ,mQ,mQ) is the mass 
matrix. The model is not renormalizable and we use a sharp cut off A to treat the divergent integrals. Following 
m we take: A = 602.3 MeV, GA^ = 1.835, KA"^ = 12.36, m^ = 5.5 MeV, mfj = 5.5 MeV, mg = 140.7 MeV. These 
parameters have been determined by fitting the /^ (pion decay constant), 711^^ (pion mass), niK (kaon mass) and m ' 
(j]' mass) to their experimental values. Due to the interaction, the mass of the quarks in the medium (dynamical 
masses) arc in general different from the bare masses (current masses). In the Hartrce approximation, the dynamical 
quark masses are determined by the solution of the gap equation: 

rrii = ml~ AG < qiqi > +2K < q^q^ X ququ >, i^ j ^k. (16) 

The quark condensates < gig,; > at zero temperature are given by: 

< m^ >= -\ r dke 7' (17) 

where kp. = (Tr^n^)^ " is the Fermi momentum of the quark i and n^ its number density. The energy density (e) and 
pressure (P) of the system are given by 
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J2 - dkk'^m^ + k^ + B,ff, (18) 



P=-e+ Y. n,^m^ + kl^. (19) 

In the NJL model the bag constant is not a phcnomenological input parameter, like in the MIT bag model. However, 
one can still define an effective bag pressure, Pe// = Po — P, generated dynamically with origin in the spontaneous 
breaking of chiral symmetry, where 
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and Bq = B{nu = n^ = rig = 0) is introduced to ensure zero pressure at zero density and temperature. 

The system of equations ([T5|) . ([TT]) is solved numerically for a fixed value of the baryonic density pB = {nu+nd+ns)/i 
and the pressure and energy density are then calculated from equations p8)) and (|T9)) . 

At finite temperature T , the above expressions can be generalized starting from the grand canonical potential per 
unit volume J7 [20|: 

f2 = ^ f2m^ + 2G < gi(7i >^ -4i^ < uu >< (Jrf >< ss > +Bo , (21) 

where 



O I 9 



The particle density of the quark i is given by: 

3 

-2 



(22) 



= ^l dk e[f,{k) - Uk)] , (23) 
while quark condensate becomes: 

< m. >=\ j dk k'^L=[Mk) + Mk) - 1], (24) 

being /j;(fc) and /i(fc) the quark and antiquark Fermi distribution functions, respectively: 

U{k) = [l + elV^'+^-^O/^') ' (25) 

/.(fc)=fl + e(^/^'+^+^0/^y' . (26) 



The gap equations are derived minimizing 17 with respect to the constituent quark masses vni. The expressions 
obtained are identical to Eqs. p^ with the quark condensate given by Eq. ([21)) . The pressure and entropy density 
of the system are given by: 



Finally the the energy density reads: 



e = -P + ^ Ui n^ + sT . (28) 



C. The Chromo Dielectric Model 

The Chromo Dielectric model |19| is a confinement model that has been extensively used to study properties of 
single nucleons or to investigate quark matter in neutron stars [35l - l43 | and supernovae explosions J43l - |45] . Confinement 
is achieved through the introduction of a scalar- isoscalar chiral singlet field x- The Lagrangian density of the model 
reads: 

where ipf represents the quark field of flavor /, t/((j, tt) is a mcxican-hat potential 
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U{a,n) = ^{a'+n^-f^), (30) 



and, for x wc consider the simplest potential 

V{x) = \mIx^ . (31) 

The characteristic feature of the CDM is that quark masses rescale as an inverse power of the field x f-^d, therefore, 
acquire a density dependence 

mu,d = 7 — : rus = — . (32) 

XU X 

In vacuum x vanishes thus providing confinement. The coupling constants are given by g^^d — gif-n i ^3) and 
Qs = g{2ffc — /tt), where /^ = 93 MeV and /if = 113 MeV are the pion and the kaon decay constants. The other two 
constants ^3 and m^ are fixed in such a way that: ^3 = /a'± — fa" = —0.75 MeV and m^ = 1.2 GeV. Following Ref. 
|43| . for the mass M^ of the field x and its coupling g we take M^ = 1.7 GeV and g = 23 MeV, respectively. These 
values lead to reasonable values for the average delta- nucleon mass and for the nucleon isoscalar radius [43 ■ 
At zero temperature the energy density in the CDM is given by: 

' = ^ E / "' '^'^ k'^Jk'+mU<y,x)+ B^ff, (33) 
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with 

Beff = V{x) + Uia,7f^0). (34) 

In mean field approximation, the field 7? vanishes while the other two fields a and x ^tre replaced by their mean values. 
The equations of motion arc obtained by minimizing the energy density of the system. One gets: 
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The inclusion of temperature in the CDM has been carried out following Ref. [36| . We can start from the grand 
canonical potential per unit volume il : 
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where rrii and fit are the effective mass (defined in Eq. ([32])) and the chemical potential of the quark i. Minimizing 
f2 with respect to a and X; we get the correspondig equations of motion at T ^ 0. 

dn\ fdn\ , , 

^ =0, — =0. 39 

The particle densities of the three quarks are given by Eq. (1231) and the pressure and energy density can be calculated 
using Eqs. ^ and ^. 



III. PHASE EQUILIBRIUM 

In the region of high density (high baryon chemical potential) and low temperature (which is the one relevant 
for neutron st ar p hysics) many QCD-inspired models suggest the deconfinement transition to be a first-order phase 
transition [4^, |40|)- Under this assumpion the conditions for phase equilibrium are thus given by the Gibbs phase rule 

Th = Tq = T, Ph = Pq^Po (40) 

^iHiT,Po)^^^QiT,Po) (41) 

where 

Eh + Ph - srT Eq+Pq-sqT 
Mff = , MQ = — ^ 42) 

are the Gibbs energies per baryon (average chemical potentials) for the hadron (H) and quark (Q) phase respectively, 
£h (eq), Ph {Pq), sh (sq) and nn {nq) denote respectively the total {i.e., including leptonic contributions) energy 
density, total pressure, total entropy density, and baryon number density of the two phases. The pressure Pq defines 
the transition point. For pressures above Pq the hadronic phase is metastablc, and the stable quark phase will appear 
as a result of a nucleation process. 

Small localized fluctuations in the state variables of the metastable hadronic phase will give rise to virtual drops of 
the stable quark phase. These fluctuation are characterized by a time scale Vq ~ 10~^^ s set by the strong interaction 
that is responsible for the deconfinement phase transition. This time scale is many orders of magnitude smaller than 
the typical time scale set by weak interactions, therefore, quark flavor must be conserved during the deconfinement 
transition. We will refer to this form of dcconfined quark matter, in which the flavor content is equal to that of the 
/3-stable hadronic system at the same pressure and temperature, as the Q* phase. Soon afterward a critical size drop 
of quark matter is formed, the weak interactions will have enough time to act, changing the quark flavor fraction of 
the deconfined droplet to lower its energy, and a droplet of /3-stable quark matter is formed (hereafter the Q phase). 

This first seed of quark matter will trigger the conversion [50l452| of the pure hadronic star to a quark star (hybrid 
star or strange star). Thus, pure hadronic stars with values of the central pressure larger than Pq are metastablc to the 
decay (conversion) to hybrid stars or to strange stars [i|tM| • The mean lifetime of the metastable stellar configuration 
is related to the time needed to nucleate the first drop of quark matter in the stellar center and depends dramatically 
on the value of the stellar central pressure [l|-lg and central temperature [lj> [l3 ■ 

In order to explore the astrophysical implications of quark matter nucleation, following Ref. [l|-|3|, we introduce 
the concept of critical mass Mcr for the hadronic star sequence. 

In the case of cold and deleptonized stars, Mcr can be defined as the value of the gravitational mass of the metastable 
hadronic star for which the nucleation time t takes a "reasonable small" value in comparison with typical ages of 
young pulsars as, e.g. the Crab pulsar. Thus, according to Ref. [I-Q we take Mcr{T = 0) = M{t = 1 yr, T = 0). It 
is worth recalling that the nucleation time is an extremely steep function of the hadronic star mass [l|-Q , therefore 
the exact value of r chosen in the definition of Mcr{T ~ 0) is not crucial. We have verified that changing r from 1 yr 
to 10^ s modifies Mcr{T = 0) by - 0.02%. 

In the case of newly formed compact stars, the characteristic evolutionary time-scale is the proto-hadronic star 
cooling time tcooi, i-C- the time it takes the new born star to reach a cold and d elep tonized configuration. The cooling 
time has been evaluated to be [S^l tcooi ~ a few 10^ s. Thus, according to Ref. [3, [I3, we consider isoentropic stellar 
configurations (with an entropy per baryon S in the range 1-2 ks) and define the critical mass for proto-hadronic 
stars as Mcr{S) = Mir = 10^ s, S* = const). 

Notice that pure hadronic stars with Mhs > Mcr are very unlikely to be observed, while pure hadronic stars with 
Mhs < ^cr are safe with respect to a sudden transition to quark matter. Thus Mcr plays the role of an effective 
maximum mass for the hadronic branch of compact stars (sec discussion in Ref. [3|). While the Oppenheimer-Volkov 
maximum mass is determined by the overall stiffness of the equation of state for hadronic matter, the value of Mcr 
will depend in addition on the properties of the intermediate non /3-stable Q* phase. 

IV. QUARK MATTER NUCLEATION RATES 

The nucleation process of quark matter in hadronic stars can take place during different stages of their evolution 
[151 . This is due to the fact that nucleation can proceed both via thermal activation or quantum tunnelling (at 
zero or finite temperature). The core of a newborn neutron star reaches temperatures of 10 — 40 MeV [3ll |32| and. 



consequently, this era is dominated by the thermal nucleation regime; on the other hand a cold delcptonized neutron 
star can nucleate quark matter only via quantum tunnelling because the thermal nucleation time diverges in the limit 
of zero temperature (see [ij, [T^ and the following discussion) . 

The energy barrier, which represents the difference in the free energy of the system with and without a (5*-matter 
droplet, can be written as 

C/(7^, T) = -TTHQ, {pLQ. - ^iH)n'^ + 47ra7e2 (43) 

where TZ is the radius of the droplet (supposed to be spherical), and a is the surface tension for the surface separating 
the hadronic phase from the Q* phase. The energy barrier has a maximum at the critical radius TZc = 2a /[nq* {fj,H — 
MQ*)]- Notice that we have neglected the term associated with the curvature energy and also the terms connected 
with the electrostatic energy, since they are known to introduce only small corrections [3|, |53| . The value of the surface 
tension a for the interface separating the quark and hadronic phase is poorly known, and the values typically used in 
the literature range within 10 — 50 MeV fm~^ |53l - [5a |. In the following, we assume a to be temperature independent 
and we take cr = 30 MeV fm^^. 

The quantum nucleation time Tq (at zero and finite temperature) can be straightforwardly evaluated within a 
semi-classical approach [5j, [5a,[53- First one computes, in the Wentzel-Kramers-Brillouin (WKB) approximation, 
the ground state energy Eq and the oscillation frequency i^o of the drop in the potential well U{TZ,T). Then, the 
probability of tunnelling is given by 

AiEo) 



Po = exp 



(44) 



where A{E) is the action under the potential barrier, which in a relativistic framework reads 
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A{E) = - / V[2^«Wt^+^^^^^C^(^¥[t^(^T^^ (45) 

being Tl± the classical turning points and m{TZ) the droplet effective mass. The quantum nucleation time is then 
equal to 

Tg - {i^oPoNcr^ , (46) 

with Nc '^ 10^® being the number of nucleation centers expected in the innermost part (r < Rnuc ^ 100 m) of the 
hadronic star, where pressure and temperature can be considered constant and equal to their central values. 
The thermal nucleation rate can be written |57| as 

/ = ^^0 exp(-[/(7ee, T)/T) (47) 

where the statistical prefactor (see Ref. [S^), is given by: 

S,Q is the quark correlation length, which gives a measure of the thickness of the interface layer between the two phases 
(the droplet "surface thickness"). In the present calculation we take ^g = 0.7 fm according to the estimate given in 
Refs. [5i,[5|. 

For the dynamical prefactor we have used a general expression which has been derived by Venugopalan and Vischer 
[HI (see also Refs. ^E^) 
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7^3(Au;)2 



XT + 2(^f^ + C 



(49) 



where Au; = wq^, — wh is the difference between the enthalpy density of the two phases, A is the thermal conductivity. 



and r] and C are, respectively, the shear and bulk viscosities of hadronic matter. According to the results of Ref. j62| . 
the dominant contribution to the prefactor k comes from the shear viscosity rj. Therefore, we take A and C equal to 
zero, and we use for the shear viscosity the following relation [62|: 

_ 7.6 X 1026 /nff\2 MeV 
"^ " (r/MeV)2 [^) fa^' ^ > 

with riQ = 0.16 fm"'^ being the saturation density of normal nuclear matter. The thermal nucleation time rj/i, relative 
to the innermost stellar region {Vnuc = i^^ / ^) ^riuc) where almost constant pressure and temperature occur, can thus 
be written as nh = (Kuc^)"^- 



V. RESULTS AND DISCUSSION 

In Fig. [l]we show the Gibbs energy per baryon for the /3-stable hadronic phase (continuous hnes) and the Q*-pha.se 
(dashed hnes) as a function of the pressure for various parametrizations of the hadronic equation of state and the two 
models, NJL (left panel) and CDM (right panel), used to describe the deconfined phase. Above the phase equilibrium 
pressure Pq the /3-stable hadron phase is metastable, and the formation of the stable (with respect to the strong 
interactions) Q*-phase will occur via a nucleation process. An interesting difference exists between the NJL and the 
Chromo Dielectric model: starting with a hadronic phase made of /3-stable nucleonic matter (continuous line labeled 
np) and next including hyperons with an increasing concentration (continuous lines labeled respectively Xa = 0.7, 
and Xa = 0.6), the transition point Pq moves to lower values in the CDM, whereas in the NJL model the opposite 
behaviour is observed. This ultimately can be traced back to the different numerical values and density (pressure) 
dependence of the dynamical strange quark mass in the two models. To elucidate this connection, we plot in Fig.[2]the 
masses (upper panels) and chemical potentials (lower panels) of the u, d and s quarks in the non-/3-stable (5*-phase 
as function of the pressure for both models and for two different parametrizations of the hadronic phase {xa- = 0.6, 
and 0.7). Here we note that the strange quark mass and chemical potential arc much larger in the NJL model than 
in the CDM one in all range of the pressure explored, and particularly at the phase transition pressure Pq. 

It is important to note that a) at P=0 the density of quarks u and d is not zero, pd ^ 2pu. Therefore the chemical 
potential is larger than the effective mass and it is larger for the d quarks; b) for P < 20-30 MeV fm~^ the density of 
s-quarks is zero and the s quark chemical potential coincides with the effective s quark mass, the chemical potential 
decreases with the pressure because the effective mass decreases; c) above P > 20-30 MeV fm~'^ the s quark density 
is nonzero and the chemical potential increases with the pressure. The s-quark chemical potential is smaller than the 
u and d quark chemical potential, although the mass is larger, because the u and d quark densities are much larger. 
At 350 MeV fm"'^ the fractions of quarks u, d and s are respectively 0.33, 0.45 and 0.22. Comparing with the NJL 
model, the big difference is the much larger effective mass the s-quark has in NJL: the onset of the s-quark occurs at 
the same pressure (set by the hyperon threshold in the /3-stable hadronic phase) but the mass is much larger. 

In order to explain the different behaviour of the two chiral approaches, we observe that there are two opposite 
effects that define the phase transition: a smaller x^ gives rise to a larger strangeness fraction in the hadronic phase 
and, therefore, to a softer hadronic EOS, shifting possible deconfincmcnt transition point to larger densities in /3-stable 
matter. On the other hand, the Q* phase has the same strangeness content of the hadronic one: in the CDM model 
all quark masses are similar and a more symmetric uds quark matter is energetically favoured, while in the NJL model 
the s quark mass is much larger and a larger ud quark fraction is favoured. This is confirmed comparing the quark 
chemical potential in both models. One can see that in the NJL model an increase of the hyperonic content of the 
hadronic phase gives rise just to a small decrease of the strange quark chemical potential, much smaller than the one 
occuring in the CDM model. 

It is also useful to analyse the effective bag pressure -Be// for NJL and CDM models as a function of the pressure 
(Fig. [3]). The value of i?e// at the transition pressure Pq is denoted with a filled circle. Coherently with the previous 
discussion, it is seen that for the NJL model the larger the hyperon content in the /3-stable hadronic phase the larger is 
the value of Pe// at the transition point. The opposite behaviour occurs for the CDM model. This is easly understood 
considering what happens in the MIT bag model for different values of the bag constant P. In this case Pe// is a 
constant and the discussion carried out in this framework is more transparent. In Fig. H] the Gibbs energy per baryon 
for Q* (thin lines) and hadronic (thick lines) matter is plotted for Xo- = 0.6 and 0.7 and three values of the bag 
pressure P. Some conclusions are immediate: i) the larger the bag value the larger the Gibbs energy per baryon for 
the Q* phase therefore the transition is disfavoured, ii) the smaller the Xa value the smaller is the Gibbs energy per 
baryon both in the hadronic and in the Q* phase, iii) the last effect is much larger in the Q* phase for the smaller 
pressures, but becomes of the same order of magnitude in both phases for P > 150 MeV fm~'^. Taking together effect 
i) and iii) and observing the very different value of BeffiPo) in the two chiral approaches considered, we deduce that 
in the NJL model the phase transition is disfavoured for the smaller value of x^ while the opposite occurs for the 
CDM model. 

The phase equilibrium curve Po{T) between the hadron and Q* phases is shown in Fig. [5] in the case of the NJL 
model (left panel) or the CDM (right panel) to describe the deconfined phase. For the hadron phase we take Xo- = 0.7 
in both cases. The region of the Pq-T plane above each curve represents the deconfined Q*-phase. We have checked 
that in the CDM, the phase separation line moves downwards (upwards) if an hadronic EOS with a larger (smaller) 
strangeness content is used, e.g., x^ = 0.6 {np matter). Using the NJL model the opposite occurs due to the reasons 
discussed above for the case of zero temperature phase transition. 

Now we can proceed to analyse the effects of quark matter nucleation process in the core of neutron stars. We 
first consider quark matter nucleation in cold /3-stable hadronic matter that is the typical situation in the core of a 
neutron star a few minutes after its birth. At this point of our discussion, we assume that quark matter nucleation 
has not occured during the proto- neutron star stage (see below). 
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In Fig. |6]we show the gravitational mass versus central pressure for different compact star models. Hadronic star 
sequences are calculated using the GMl parametrization for pure nucleonic matter (black curve), hyperonic matter 
with Xa = 0.7 (blue curve) and x^ = 0.6 (red curve). The quark star (QS) sequence is represented by the green curve. 
For the quark models cosidered in this paper, all QS sequences are made of hybrid stars (YS). Results in the left 
(right) panel are relative to the NJL (CDM) model for the quark phase. The configuration marked with an asterisk 
represents, in all cases, the hadronic star for which the central pressure is equal to Pq and thus the quark matter 
nucleation time is r = oo. The critical mass configuration is denoted by a full circle. The stellar conversion process 
llri^] of the critical mass configuration into a final quark star with the same stellar baryonic mass (filled square) is 
denoted by the dashed line connecting the the circle to the square. Notice that in most of the cases reported in the 
figure the quark matter nucleation process will lead to the formation of a black hole (for these cases we do not plot in 
Fig.inithe coresponding YS sequence). In particular, within the present values for the EOS parameters the formation 
of quark stars is not possible modeling the quark phase with the NJL model. 

We next consider the case of new born hadronic stars (proto-hadronic stars, PHSs). In this case the quark decon- 
finement phase transition is likel y tr iggered by a thermal nucleation process and it will occur in those PHSs with a 
gravitational mass M > Mcr{S) [IJ, [la]- Here we consider the case of neutrino free matter, since it has been shown 
|15| that neutrino trapping does not alter substantially the outcomes of the PHS evolution. 

The evolution of a PHS within this scenario is delineated in Fig. [71 where we plot the appropriate stellar equilibrium 
sequences in the gravitational-baryonic mass plane obtained from the CDM for the quark phase and the GMl model 
in the case of hyperonic matter with x^ = 0.6 (left panel), and x^ — 0.7 (right panel). In particular, we plot the PHS 
sequence, i.e. isoentropic hadronic stars {S = 2 ks) and neutrino- free matter (upper line), and the cold hadronic 
star (HS) sequence (dashed line). The asterisk and the full circle on these lines identify respectively the stellar 
configuration with r = oo and the critical mass configuration. We denote as Mg^f = MB,cr{S = 2kB) the baryonic 
critical mass for the PHS sequence and as Mg^^ = MB,cr{S ~ 0) the baryonic critical mass for the cold hadronic 
star sequence. The lower continuous (green) line represents the cold QS sequence having a a maximum gravitational 
(baryonic) mass M^f^ i^B max)- ^^ assume [53 ( Mb = constant during these stages of the steUar evolution ). We 
note, indeed, that sizeable mass accretion on the proto-neutron star occurs within a time of ~ 0.5 s after core bounce 
[3ll[32])- During the subsequent stages, the star thus evolves with Mb — const. . 

Thus according to the results in the left panel of Fig. [71 proto-hadronic stars with a baryonic mass AIb < M^^J' 
(— 1.16 il/0 within the selected EOS parametrization) will survive Q*-matter early nucleation (i.e. nucleation within 
the cooling time tcooi ~ a few 10^ s) and in the end will form stable (r — oo) cold hadronic stars. Proto-hadronic stars 
with M§^f < Mb < M^^^^^ (= 1.79 Mq for the present EOS) wiU experience early nucleation of a Q*-matter drop 

and will ultimately form a cold deleptonized quark star. The last possibility is for PHSs having Mb > M^^^^. In 
this case the early nucleation of a Q*-matter drop will trigger a stellar conversion process to a cold QS configuation 
with Mb > Mg ^^^. thus these PHSs will finally form black holes. A similar evolutionary path is found in the case 
of Xa = 0.7 (right panel). 

In Fig. [SI we plot the PHS, cold HS, and cold YS sequences in the gravitational-baryonic mass plane for the case of 
the NJL model for the quark phase and the GMl model in the case of pure nucleonic matter (right panel) or hyperonic 
matter with Xa = 0.7 (left panel). It is clearly seen that in the case of the NJL model it is almost impossible to 
populate the YS branch. Cold quark stars can be formed in the case of x^ ~ 0.7 (left panel) for a very narrow range 
of baryonic stellar masses 2.20 < Mb/Mq < 2.23M0. 

Finally, in tables [J and [IT] we report the values of the gravitational and baryonic critical mass for the PHS and HS 
sequences for the two adoped quark matter models. As it has been found in previous works [IJ, [lalj thermal effects 
reduce the values of the critical mass and increase the portion of the quark star branch that can populated via the 
stellar conversion process [l|, [2, |52| . Notice that the maximum mass configuration of the YS sequence is insensitive, 
in the case of the CDM, to the value of the hyperon coupling Xa, since in this case the threshold density for quark 
deconfinement is much lower than the density for the onset of hyperons. 



VI. CONCLUSION 

In this paper we have studied the nucleation of quark matter in both hot and cold /3-stable hadronic matter using two 
different models with chiral symmetry to describe the quark phase: the Nambu-Jona-Lasinio model and the Chromo 
Dielectric model. For the hadronic phase we chose the GMl parametrization of the non-linear Walecka model and 
we have considered pure nucleonic matter as well as hyperonic matter with a large hyperon fraction (xcr=0.6), and a 
small hyperon fraction (xo- = 0.7). 

The nucleation process forms a short-lived transitory phase (Q* phase) which has the same flavor content of the 
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initial /3-stablc hadronic phase. This particular circumstance, together with the different pressure (density) dependence 
of the strange quark effective mass in the two employed quark matter models produces considerable differences on 
the bulk properties of the phase transition and on neutron star composition and early evolution. More precisely, we 
found that for the NJL model the presence of hyperons disfavour the phase transition pushing the transition point to 
very high densities while with the CDM the opposite behaviour has been observed. In addition, we found that in the 
case of the NJL model it is almost impossible to populate the quark star branch and that quark matter nucleation 
will lead to the formation of a black hole. Thus within the NJL model for the quark phase, all compact stars are pure 
hadronic stars. In the case of the CDM, thermal effects reduce the value of the critical mass, and both hadronic and 
quark star configurations can be formed as a result of the evolution of proto-hadronic stars, depending on the value 
of the stellar baryonic mass. 

A very recent measurement [63 of the mass of the pulsar PSR J1614-2230 makes it the most massive neutron star 
known to date with a mass M = (1.97 ± 0.04) Mq. Within the EOS models employed in the present work, the 
compact star in PSR J1614-2230 could only be a pure HS (in the case the quark phase is described by the NJL model. 
Fig. [S]) formed from the evolution of a PHS with initial baryonic mass Mb < M^^f and after a long-term mass 
accretion {Maccr ~ 0.1 - 0.2 Mq) from a companion star in a binary system. This long term evolution can finally 
form pure hadronic star with a mass M < M^/ , with M^^-^ = 2.025AfQ (case x„ = 0.7) or M^/ = 2.238M0 (case 
of np matter). The CDM model fails to predict such a high mass. 
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M^/ 


A/f.i 


aC"^ 


t.fPHS 




"'J-B,max 


x^ = 


= 0.7 


2.025 


2.342 


1.964 


2.201 


1.943 


2.293 


np 


2.238 


2.634 


2.112 


2.386 


1.988 


2.287 



TABLE I. Gravitational M^ {M^r^ ) and baryonic M^f^ (M^^f) critical mass for the cold hadronic star (proto-hadronic 
star) sequence. The gravitational and baryonic maximum mass for the cold hybrid star sequence are denoted respectively as 
^max s-nd Mg^^^. The values of stellar masses are in unit of the solar mass (M© — 1.989 x 10"^"^ g). All the results reported 
in the table are relative to the GMl equation of state for the hadronic phase and the Nambu-Jona-Lasinio (NJL) model for 
the quark phase. The gravitational maximum mass for the cold hadronic star sequence is M,^^^ — 2.O42M0 in the case of 
hyperonic matter with x^ = 0.7, and M^^x ~ 2.364M© in the case of nucleonic (np) matter. 
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1.557 


1.793 



TABLE IL Same as in table[T]but using the CDM for the quark phase. The gravitational maximum mass for the cold hadronic 
star sequence is M^^x ~ 1.791A/0 in the case x^ = 0.6, and M^ax ~ 2.O42M0 in the case Xa = 0.7. 
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FIG. 1. (Color online) Gibbs energy per baryon at zero temperature (T = 0) as function of the pressure for the hadronic 
(continuous lines) and the Q* (dashed lines) phase. The arrows and the corresponding numbers indicate the value of the 
pressure Po at which the bulk phase transition takes place. Results for the NJL (CDM) model are shown in the left (right) 
panel. See text for more details. 
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FIG. 2. (Color online) Quark masses (upper panels) and chemical potentials (lower panels) in the Q* phase at T = using 
the NJL (left panels) and the CDM (right panels) models to describe the quark deconfined phase. Results for Xa = 0.6 and 
Xa — 0.7 are shown by the red continuous lines and the blue dashed ones, respectively. 
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FIG. 3. (Color online) Effective bag pressure at T = for the NJL (left panel) and the CDM (right panel) models. Results for 
Xcr ~ 0.6, Xcr = 0.7 and the pure nucleonic matter case (np) are shown. The filled circle denotes the value of the effective bag 
pressure at the transition point Pq. 
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FIG. 4. (Color online) Gibbs energy per baryon at T = for the /3-stable hadronic phase, with Xa = 0.6 (continuous thick line) 
and Xa ~ 0.7 (dashed thick line), and for the Q* phase (thin lines). The MIT bag model has been used for the Q* phase with 
a bag constant B = 80, 100 and 165 MeV/fm^. 
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FIG. 5. (Color online) Phase equilibrium curve between the /3-stable hadronic phase, with lo- ~ 0.7, and the Q* phase using 
the NJL (left panel) and the CDM (right panel) models. 
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FIG. 6. (Color online) Gravitational mass versus central pressure for compact stars. Hadronic star sequences are calculated 
using the GMl parametrization for pure nucleonic matter (black curve), hyperonic matter with Xa = 0.7 (blue curve) and 
Xa = 0.6 (red curve). The hybrid star (YS) sequence is represented by the green curve. Results in the left (right) panel are 
relative to the NJL (CDM) model for the quark phase. The configuration marked with an asterisk represents in all cases the 
hadronic star for which the central pressure is equal to Po and thus the quark matter nucleation time is r = cxd. The critical 
mass configuration is denoted by a full circle. The conversion process of the critical mass configuration into a final quark star 
with the same stellar baryonic mass (filled square) is denoted by the dashed line connecting the the circle to the square. In 
most of the cases reported in the figure the quark matter nucleation process will lead to the formation of a black hole (go to 
BH). For these cases we do not plot the coresponding YS sequence. 
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FIG. 7. Evolution of a proto-hadronic star in the gravitational-baryonic mass plane using the Chromo Dielectric model (CDM) 
to describe the quark phase and the GMl model with x^ ~ 0.6 (left panel) and Xcr = 0.7 (right panel) for the hadronic phase. 
The upper (red) line represents the stellar equilibrium sequence for neutrino-free proto-hadronic stars (PHS) with S = 2 ks- 
The (black) dashed line represents the cold the HS sequence. The asterisk and the full circle on these lines represent the stellar 
configuration with nucleation time r = oo and the critical mass configuration, respectively. The lower (green) line represent 
the cold YS sequence. Assuming Afs = const, the evolution of a PHS in this plane occurs along a vertical line. 
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FIG. 8. Same as in the previous figure but in the case of the NJL model for the quark phase and the GMl model for hyperonic 
matter with Xa =0.7 (left panel) and pure nucleonic (np) matter (right panel). 



